Three-dimensional finite element simulation of fracture propagation in rock specimens with pre-existing fissure(s) under compression and their strength analysis     

Ramin Pakzad  Shanyong Wang  Scott W. Sloan 《国际地质力学数值与分析法杂志》2020,44(10):1472-1494

A FORTRAN program, consistent with the commercially available finite element (FE) code ABAQUS, is developed based on a three-dimensional (3D) linear elastic brittle damage constitutive model with two damage criteria. To consider the heterogeneity of rock, the developed FORTRAN program is used to set the stiffness and strength properties of each element of the FE model following a Weibull distribution function. The reliability of the program is assessed against available experimental results for granite cylindrical specimens with a throughgoing, flat and inclined fissure. The calibration procedure of the material parameters is explained in detail, and it is shown that the compressive to tensile strength ratio can have a substantial influence on the failure response of the specimens. Numerical simulations are conducted for models with different levels of heterogeneity. The results show a smaller load bearing capacity for models with less homogeneity, representing gradual coalescence of fully damaged elements forming throughout the models during loading. The maximum load bearing capacity is studied for various combinations of inclination angles of two centrally aligned, throughgoing and flat fissures of equal length embedded in cylindrical models under uniaxial and multiaxial loading conditions. The key role of the compressive to tensile strength ratio is highlighted by repeating certain simulations with a lower compressive to tensile strength ratio. It is proven that the peak loads of the rock models with sufficiently small compressive to tensile strength ratios containing two throughgoing fissures of equal length are similar, provided that the minimum inclination angles of the models are the same. The results are presented and discussed with respect to the existing experimental findings in the literature, suggesting that the numerical model applied in this study can provide useful insight into the failure behaviour of rock-like materials.  相似文献   

Effects of the geometry of two‐dimensional fractures on their hydraulic aperture and on the validity of the local cubic law          下载免费PDF全文

Ahmad Hajjar  Luc Scholtès  Constantin Oltéan  Michel Antoine Buès 《水文研究》2018,32(16):2510-2525

Flow through rough fractures is investigated numerically in order to assess the validity of the local cubic law for different fracture geometries. Two‐dimensional channels with sinusoidal walls having different geometrical properties defined by the aperture, the amplitude, and the wavelength of the walls' corrugations, the corrugations asymmetry, and the phase shift between the two walls are considered to represent different fracture geometries. First, it is analytically shown that the hydraulic aperture clearly deviates from the mean aperture when the walls' roughness, the phase shift, and/or the asymmetry between the fracture walls are relatively high. The continuity and the Navier–Stokes equations are then solved by means of the finite element method and the numerical solutions compared to the theoretical predictions of the local cubic law. Reynolds numbers ranging from 0.066 to 66.66 are investigated so as to focus more particularly on the effect of flow inertial effects on the validity of the local cubic law. For low Reynolds number, typically less than 15, the local cubic law properly describes the fracture flow, especially when the fracture walls have small corrugation amplitudes. For Reynolds numbers higher than 15, the local cubic law is valid under the conditions that the fracture presents a low aspect ratio, small corrugation amplitudes, and a moderate phase lag between its walls.  相似文献   

Solution for a plane strain rough‐walled hydraulic fracture driven by turbulent fluid through impermeable rock          下载免费PDF全文

N. Zolfaghari  E. Dontsov  A. P. Bunger 《国际地质力学数值与分析法杂志》2018,42(4):587-617

The impact of turbulent flow on plane strain fluid‐driven crack propagation is an important but still poorly understood consideration in hydraulic fracture modeling. The changes that hydraulic fracturing has experienced over the past decade, especially in the area of fracturing fluids, have played a major role in the transition of the typical fluid regime from laminar to turbulent flow. Motivated by the increasing preponderance of high‐rate, water‐driven hydraulic fractures with high Reynolds number, we present a semianalytical solution for the propagation of a plane strain hydraulic fracture driven by a turbulent fluid in an impermeable formation. The formulation uses a power law relationship between the Darcy‐Weisbach friction factor and the scale of the fracture roughness, where one specific manifestation of this generalized friction factor is the classical Gauckler‐Manning‐Strickler approximation for turbulent flow in a rough‐walled channel. Conservation of mass, elasticity, and crack propagation are also solved simultaneously. We obtain a semianalytical solution using an orthogonal polynomial series. An approximate closed‐form solution is enabled by a choice of orthogonal polynomials embedding the near‐tip asymptotic behavior and thus giving very rapid convergence; a precise solution is obtained with 2 terms of the series. By comparison with numerical simulations, we show that the transition region between the laminar and turbulent regimes can be relatively small so that full solutions can often be well approximated by either a fully laminar or fully turbulent solution.  相似文献   

Micromechanical approach to effective viscoelastic properties of micro‐fractured geomaterials          下载免费PDF全文

Cássio B. Aguiar  Samir Maghous 《国际地质力学数值与分析法杂志》2018,42(16):2018-2046

The aim of this paper is to formulate a micromechanics‐based approach to non‐aging viscoelastic behavior of materials with randomly distributed micro‐fractures. Unlike cracks, fractures are discontinuities that are able to transfer stresses and can therefore be regarded from a mechanical viewpoint as interfaces endowed with a specific behavior under normal and shear loading. Making use of the elastic‐viscoelastic correspondence principle together with a Mori‐Tanka homogenization scheme, the effective viscoelastic behavior is assessed from properties of the material constituents and damage parameters related to density and size of fractures. It is notably shown that the homogenized behavior thus formulated can be described in most cases by means of a generalized Maxwell rheological model. For practical implementation in structural analyses, an approximate model for the isotropic homogenized fractured medium is formulated within the class of Burger models. Although the approximation is basically developed for short‐term and long‐term behaviors, numerical applications indicate that the approximate Burger model accurately reproduce the homogenized viscoelastic behavior also in the transient conditions.  相似文献   

浙西北江山─绍兴断裂带构造演化特征     

竺国强  姜继双  陈梓军  周培华 《浙江国土资源》1997,(2)

江山—绍兴断裂带是中国南方最重要的构造带之一。它是江南陆块与华夏陆块神功运动形成的碰撞拼贴带。它表现出明显的多期变形、构造叠加的特征：神功运动形成近东西向平卧褶皱；印支运动沿断裂带有向北西强烈的推覆挤压。断裂带不同时期的活动反映了不同特征的构造环境，形成不同的构造变形特征，并对两侧的地质历史具有明显的控制作用。  相似文献   

The Geological-Geophysical-Geochemical Exploration Model of the Nancha Type Gold Deposit in Tonghua Area, Jilin Province     

《东北亚地学研究》2001,(2)

  相似文献   

A crack model of afterslip events on shallow faults     

Maria Elina Belardinelli  Maurizio Bonafede 《Geophysical Journal International》1991,106(3):521-530

  相似文献   

元岭金矿蚀变破碎带特征及其找矿意义     

刘红樱  周顺之  胡受奚  王群昌 《河南国土资源》1996,(2)

根据对元岭金矿三条蚀变破碎带一般地质特征、微量元素特征的对比研究，并借鉴邻近矿区蚀变破碎带的矿化特征，本文阐明了元岭蚀变破碎带矿化的有利构造部位，总结了蚀变破碎带特征与矿化强度之间的关系，指出蚀变破碎带的矿化蚀变元素为中低温组合，矿化蚀变愈强则参与元素的种类愈多、数量愈大。围岩性质对矿化蚀变有一定影响，围岩中的易迁移亲铜元素和亲硫铁族无素对矿化有利。据此，为进一步成矿预测、评价蚀变破碎带的含矿性提供依据。  相似文献   

Numerical modelling of the impact crater depth-diameter dependence in an acoustically fluidized target     

K. Wünnemann  B.A. Ivanov 《Planetary and Space Science》2003,51(13):831-845

2D numerical modelling of impact cratering has been utilized to quantify an important depth-diameter relationship for different crater morphologies, simple and complex. It is generally accepted that the final crater shape is the result of a gravity-driven collapse of the transient crater, which is formed immediately after the impact. Numerical models allow a quantification of the formation of simple craters, which are bowl-shaped depressions with a lens of rock debris inside, and complex craters, which are characterized by a structural uplift. The computation of the cratering process starts with the first contact of the impactor and the planetary surface and ends with the morphology of the final crater. Using different rheological models for the sub-crater rocks, we quantify the influence on crater mechanics. To explain the formation of complex craters in accordance to the threshold diameter between simple and complex craters, we utilize the Acoustic Fluidization model. We carried out a series of simulations over a broad parameter range with the goal to fit the observed depth/diameter relationships as well as the observed threshold diameters on the Moon, Earth and Venus.  相似文献   

基于统一强度理论的地基临界荷载公式   总被引：1，自引：0，他引：1  

范文  林永亮  秦玉虎 《地球科学与环境学报》2003,25(3):48-51

基于统一强度理论,按照临界荷载公式的推导方法,得到地基临界荷载的统一解公式。Mohr-Coulomb强度理论计算的结果为其特例。把计算结果与《建筑地基基础设计规范》中的结果进行了对比分析,认为当不同程度地考虑中间主应力σ2的影响时,可以有效地发挥地基的强度潜能。并通过实例计算了统一强度理论参数b取不同值时地基承载力特征值,并进行了对比分析。  相似文献